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EVOLUTION  OF  A  FINITE  PULSE  OF  RADIATION  IN  A 
HIGH-POWER  FREE-ELECTRON  LASER 


I.  Introduction 

Hany  of  the  analyses  and  numerical  studies  of  the  free-electron  laser 

(FEL)  have  been  directed  towards  in-depth  investigation  of  individual 

aspects  of  the  underlying  mechanism.  What  emerges  is  that  the  FEL 

interaction  is  sufficiently  complex  to  be  endowed  with  a  variety  of 

interesting  physical  phenomena.^  The  purpose  of  the  present  paper  is  to 

discuss  the  results  from  a  numerical  study  in  a  sufficiently  general 

setting  wherein  several  of  the  phenomena  of  interest  are  simultaneously 

operative.  The  particular  case  examined  is  that  of  a  high-current  (-  kA), 

high-power  (~  GW},  short  wavelength  (~  pm),  tapered-wiggler  FEL.  The 

results  indicate  that  in  addition  to  enhancing  the  extraction  efficiency, 

tapering  of  the  wiggler  has  the  further  benefit  of  improving  the  quality  of 

the  output  by  suppressing  sideband  modulation  of  the  optical  field.  The 

role  of  optical  guiding  is  emphasized.  Specifically,  the  results  are  used 

to  illustrate  the  characteristics  associated  with  refractive  guiding  and 

gain  focusing.  The  gradual  transition  between  the  two  forms  of  guiding  as 

2 

the  tapering  rate  is  varied  is  discussed  in  detail. 
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II .  Mathematical  Model 

The  purpose  of  this  section  is  to  presert  the  complete  set  of 
equations  forming  the  basis  of  the  numerical  model.  The  vector  potential 
of  the  optical  field  is  given  by 

As  =  |  exp[ig  z  -  «t)]ex  +  c.c.f 

where  A  is  the  slowly-varying  envelope,  w  is  the  radian  frequency,  c  is  the 
speed  of  light  in  vacuo,  and  ex  is  the  unit  vector  along  the  x  axis.  In 

n 

the  assumed  azimuthally  symmetric  geometry  considered,  a  =  |e  |A/mc  may  be 
expanded  as 


a(r,z,„  .  Efl  *„(*,<>  L„  (2r2/rs2) 

x  exp[-(l  -  ia)  r2/rs2J,  (1) 

where  — | e |  is  the  charge  on  an  electron  of  rest-mass  m,  Ln  is  the  Laguerre 

polynomial  of  order  n,  and  r  is  the  radial  coordinate.  In  this  expansion 

rs(z, t)  is  the  spot-size,  a(z,t)  is  related  to  the  curvature  R-^  of  the 
-1  2 

optical  field  by  R  =  2c<*/<ors,  and  an(z,t),  n  =  0,  1,  2,  ...,  are  the 
expansion  coefficients.  Substituting  Eq.  (1)  into  the  well-known  parabolic 
wave  equation  and  employing  the  source-dependent  expansion  technique  of 

4 

Ref.  3,  one  obtains  the  following  set  of  equations 
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where 
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A  =  (2n  +  1  -  ia)  +  i  (2n  B_  +  B), 

n  i  k 

tor 

s 


B  =  F1/a0, 


■  Br  +  iBj.  (2f) 

In  the  formula  for  F^,  v  =  1^  [kA]/17,  where  1^  is  the  electron  beam 

current  in  kilo-Amperes,  the  summation  is  over  the  electrons  at  radial 

locations  { r ^ }  in  a  given  ponderomotive  bucket,  N  being  the  number  of 

electrons  initially  therein;  ty,  =  (to/c  +  k  )z.  ~  ttt  -s  phase  of  the  jth 

J  W  J 

electron,  located  at  z  ^ ,  relative  to  the  ponderomotive  wave.  The 
normalized  wigg'Ur  amplitude,  ay,  is  defined  in  Eq.  (5d). 

A  detailed  presentation  of  the  source-dependent  expansion  approach  is 
given  in  Ref.  3.  For  orientation,  however,  it  should  be  noted  that  in 
vacuo  one  has  the  well-known  result 


rs(z>  -  rs(0)  (l  ♦  z2/Zj2)  , 

a(  z  )  =  z  /  , 


where  rs(0)  is  the  minimum  spot  size  (at  z  =  0)  and 
7-r  -  h  rs2<°> 


is  the  Rayleigh  range.  These  results  follow  from  Eqs.  (2)  upon  neglecting 
B,  i.e.,  neglecting  thr  electron  beam,  and  assuming  that  the  wavefronts 
are  plane  at  z=0. 
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The  scalar  potential  of  the  viggler  is  represented  by 

X  =  -  (B  /k  )  cosh  (k  x)  sinh  (k  y)  cos  (k  z),  (A) 

v  y  x  y  w 

where  B  is  the  induction  and,  for  the  parabolic  pole  faces  considered 

herein,  the  wavenumbers  along  the  three  axes  are  related  via  k  =  k  = 

x  y 

k  /V2.  With  the  scalar  potential  given  by  Eq.  (4),  tha  betatron  motion  of 

V 

the  electrons  is  described  by  the  quasi-classical  (VKB)  solutions 


xfl  .  cos 


Jdt-  <.6(f) .  d. 


and 


Ia 

y  =  cos 


J  dt'  tog(t')  +  ^j, 


where  x_,  y_,  €  ,  aid  F  are  constants, 
P  P  *  y 
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a  ck 
w  w 


3  ”  2y 


is  the  betatron  frequency,  and 


a  = 
w 


e  B. 


k  me 
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(5a) 


(5b) 


(5c) 


(5d) 


is  the  normalized  wiggler  amplitude. 

2 

Finally,  the  synchrotron  motion  of  an  electron  of  energy  ymc  located 
at  (z,r)  is  described  by 


,  i  oia  f_ 
d  y  V  B  y* 

dt  4  y  ^ 


a  L 


n  n 


exp  i  ^ 
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(l-io)r2 


+  c .  c .  , 


(6a) 
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djjj 
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2 
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(6b) 


where  ip  =  ( w/c  +  k  )  z  -  wt  is  the  phase  of  the  electron  in  the 

V 

ponderomotive  wave,  and  fg  =  JQ(0  -  J^(C)  is  the  usual^  difference  of 

2  2 

Bessel  functions  with  £,  =  (a  /2)  /(I  +  a„  /2). 

W  V 
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III.  Initialization 


The  initial  electron  distribution  is  as  follows.  The  electrons  are 

uniformly  loaded  on  the  mesh  corresponding  to  the  z  axis,  and  they  are 

weighted  in  such  a  way  as  to  represent  a  parabolic  density  profile  along  z. 

For  the  transverse  coordinates  (x,v  ,y,v  )  the  electrons  are  loaded 

x  y 

randomly  in  such  a  way  that  i)  the  radial  density  profile  is  parabolic,  and 
ii)  the  transverse  velocity  profile  is  Gaussian.  The  normalized  edge 
emittance  for  such  a  distribution  may  be  defined  in  a  manner  similar  to 
that  for  a  Gaussian-Gaussian  distribution  given  in  Ref.  7: 


1/2 

ex  =  2  /6  [<x2XY2ex2>  -  <xrf3x>2]  , 


(7) 


and  similarly  for  Here  <...>  indicates  an  average  over  the  entire 

beam.  For  a  monoenergetic  electron  beam  with  an  initial  parabolic  radial 
density  distribution 


nb(r)  =  nQ  (1  -  r2/rbo2), 


(8a) 


and  Gaussian  velocity  distribution 


x  y 


<  6‘  1 

exn  f  ^  ) 

2a  2j 

expr 

V  2a 

(8b) 


the  normalized  edge  emittance  is  given  by 


(9) 


where  r,  is  the  initial  electron  beam  radius,  a  and  a  are  the  initial 
bo  x  y 

velocity  space  widths,  and  0  =  (vx/c, v^/c, v^/c) . 

From  Eq.  (6b)  the  relativistic  factor  for  a  synchronous  (i.e., 
resonant)  electron  with  no  betatron  motion  is  given  by  d^/dt  =  0: 
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\  - I  <«'<*,>  <1  *  k2>-  <10) 

Tne  form  of  tapering  employed  in  the  computations  is  obtained  simply  by 
prescribing  a  constant  rate  of  decrease  of  energy  dyr/dz  <  0  for  a 
synchronous  electron.  With  dy^/dz  equal  to  a  given  constant,  one  obtains 
aw(z)  from  Eq.  (10). 

The  parameters  for  the  computations  presented  herein  are  listed  in 
Tables  I  and  II.  The  high-current,  high-energy  electron  beam  would 
correspond  to  an  rf  linear  accelerator  similar  to  that  at  the  Boeing 

g 

Aerospace  Company  (BAC).  It  must  be  pointed  out  that,  to  reduce  run-time 
on  the  computer,  the  electron  pulse  length  used  here  is  shorter  than  that 
at  BAC. 

The  initial  intrinsic  energy  spread  is  taken  to  be  zero.  The  last  two 
terms  in  Eq.  (6b),  which  are  due  to  wiggler  gradients,  vary  with  the 
initial  conditions  for  each  electron,  but  are  approximately  constant  along 
any  individual  orbit.  This  is,  of  course,  the  principal  virtue  of  canted 
magnetic  pole  faces. ^  Hence,  there  is  no  constraint  on  the  bucket  height 
from  wiggler  gradients  as  such.  However,  there  is  an  effective  energy 
spread  due  to  emittance  which  must,  for  a  tapered  wiggler,  be  less  than  the 
bucket  height.  From  Eq.  (6b)  it  is  simple  to  show  t..at  the  shift  6yr  in 
the  synchronous  energy  yr>  defined  by  Eq.  (10),  due  to  the  initial  spread 
in  transverse  coordinates  and  velocities  is  given  by 


P-] 


ill  *  *v2)/2 


emit 


K.)  ( 


1  +  a  /2 

w 


) 


(11) 


where  e  ,  e  and  r,  are  defined  by  Eqs.  (7)  and  (8a).  To  derive  Eq.  (11) 
x  y  do 

use  has  been  made  of  the  fact  that,  for  a  matched  electron  beam,  the 
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electron  beam  envelope  equation  implies 
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(12) 


,  ,  2.  x 

kfl<x  >  =  -  , 

P  2/6y 


2  2  2  2 
and  that,  from  (8a),  <x  >  =  r^  /6.  Here,  <x  >  indicates  the  average  of  x 

taken  with  respect  to  the  weight  function  given  by  the  density  n^(r),  Eq. 

(8a).  In  Eq.  (12)  ka,  the  betatron  wave  number,  is  related  to  the  betatron 

P 

frequency  in  Eq.  (5c)  via 


“p  =  3z<*p- 


(13) 


There  is,  of  course,  a  condition  identical  to  Eq.  (12)  relating  <y  >  to  e^. 

Retaining  the  lowest  order  optical  mode  aQ  in  Eq.  (6a)  and  dropping 
the  betatron  terms  in  Eq.  (6b),  the  full  ponderomotive  bucket  height  for 
electrons  on  axis  is  given  by 

£r]  ■  (?rr?7ri  tcos  "  8  sg"t*  _  0sin^]}  '<u> 

bucket  a  +  w'^J 


where 


ao  =  laDl  exP  (i*)* 


Sr  = 


(15) 

(16) 


Yr  is  defined  by  Eq.  (10),  and  the  phase  ^  is  related  to  the  rate  of 
change  of  rr  in  the  tapered  FEL  by 


dYr 

dt 


a  a  ufD 
w  1  o  '  B 


2r. 


sin  («(;  +  $). 


(17) 


Table  II  lists  the  full  bucket  height  for  the  various  tapering  rates 
dy^/dz  =  c'1  dyr/dt  employed  in  the  computations.  Comparing  the  set  of 
bucket  heights  with  the  effective  energy  spread  due  to  emittance, 


■8 


[ &y  /y  ]  •  ,  one  observes  that  tor  the  input  power  of  450  MU  essentially 

all  the  electrons  are  i  litially  trapped.  Also  listed  in  Table  II  are  the 
initial  resonance  phases  Eq.  (16),  for  the  various  tapering  rates. 
Note  tnat  since  in  the  computations  dy^dz  is  prescribed  to  have  a  fixed 
value  throughout  a  run,  in  general  Er  is  a  variable  through  the  viggler. 
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IV.  Results  and  Discussion 

In  this  section  numerical  results  from  the  solution  of  Eqs.  (2),  (5) 
and  (6)  are  presented.  As  mentioned  in  the  Introduction,  the  results  vill 

be  used  to  illustrate  some  of  the  key  physical  phenomena  of  the  PEL 

.  .  2  _.  10-18  .  ...  3,10,19-28 

mechanism.  These  include  sideband  growth,  optical  guiding, 

10  19  20 

efficiency  enhancement,  ’  ’  and  pulse  slippage. 

As  indicated  in  Table  II  the  conclusions  presented  in  this  paper  are 
based  on  numerical  results  that  cover  a  wide  range  of  tapering  rates.  For 
brevity,  the  results  for  only  two  tapering  rates  vill  be  discussed.  Case 
(a),  -dy  /dz  =  0.1  m-*  has  a  slow  taper,  and  Case  (b),  -dy^/dz  =  1.3 
has  a  fast  taper.  The  two  cases  discussed  in  detail  suffice  since  there  is 
a  gradual  change  in  the  physical  characteristics  in  going  from  one  tapering 
rate  to  the  next.  The  corresponding  tapering  of  the  wiggler  amplitude,  av, 

is  shown  in  Figs.  1  (a)  and  (b),  where  ay  barely  changes  for  Case  (a).  As 

shown  in  Figs.  2  (a)  and  (b),  the  resonant  angle  for  Case  (a)  starts  out  to 
be  =  2.6°  and  soon  settles  down  to  approximately  £r  -  1.5°,  while  for 
Case  (b)  the  resonant  angle  changes  from  an  initial  value  of  36.5°  to  a 
final  value  of  -  23°.  In  general,  the  resonance  phase,  e  +  Eq. 
(16),  varies  along  the  wiggler.  For  all  the  tapering  rates  indicated  in 
Table  II,  it  is  found  that  after  an  initial  transient,  £>r  settles  down  to  a 
fairly  constant  value  through  the  rest  of  the  run.  The  gain  in  radiation 
power  along  the  wiggler  is  shown  in  Figs.  3  (a)  and  (b).  The  tapering  rate 

in  Case  (b)  is  close  to  the  optimal  value  and  Fig.  3  (b)  shows  substantial 

increase  in  the  output  radiation  power. 

(a)  Sideband  Instability  and  Pulse  Modulation 

Ordinarily,  tapering  of  the  wiggler  field  is  a  means  of  enhancing  the 
extraction  efficiency  of  a  FEL.  In  previous  work  it  has  been  shown  that 
tapering  has  the  additional  benefit  of  reducing  sideband 
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modulation.  ’  ’  ’  In  particular,  it  is  shown  in  Ref.  4  that  as  the 

tapering  rate  is  increased,  the  modulation  of  the  output  signal  due  to 
sideband  frequencies  is  reduced.  Of  course,  there  is  an  optimal  rate  of 
taper,  beyond  which  excessive  electron  detrapping  from  the  ponderomotive 
buckets  leads  to  a  reduction  in  the  extraction  efficiency. 

For  the  finite  pulse  simulations,  it  also  turns  out  that  the  tapering 
rate  has  a  dramatic  effect  on  the  profile  of  the  output  optical  pulse. 

This  is  illustrated  rather  veil  in  Figs.  4  (a)  and  (b)  which  display  the 
normalized  radiation  field  amplitude  a  (z)  at  the  wiggler  exit,  and  Figs.  5 
(a)  and  (b)  which  show  the  longitudinal  profiles  of  the  power  in  the 
radiation  pulses  for  Cases  (a)  and  (b).  Figures  6  (a)  and  (b)  show  the 
transverse  profiles  of  the  radiation  pulse  for  Cases  (a)  and  (b).  Note 
that  in  all  the  profiles  for  Case  (a),  there  is  significant  amount  of 
modulation.  The  respective  spatial  Fourier  spectra  |a(X)|  as  a  function  of 
wavelength  X  are  shewn  in  Figs.  7  (a)  and  (b).  The  Fourier  spectrum  for 
the  slow  taper,  Fig.  7(a),  indicates  two  prominent  sidebands  on  either  side 
of  the  main  component,  the  anti-Stokes  sideband  in  the  vicinity  of  0.9775 
ym  and  the  Stokes  sideband  in  the  vicinity  of  0.9915  ym.  The  two  sidebands 
are  seen  to  be  of  approximately  equal  amplitude  and  of  equal  separation 
relative  to  the  carrier.  For  the  fast  taper,  Fig.  7(b),  we  notice  that  the 
sidebands  are  significantly  reduced  compared  to  those  in  Fig.  7(a).  It  is 
interesting  to  note  that  the  carrier  components  and  the  anti-Stokes 
components  have  merged  into  a  broad-band  feature  in  Fig.  7(b). 

As  is  vell-knovn,  the  wavelengths  of  the  upper  and  the  lower  sidebands 

,  .  10-14 

are  given  by 


1 1 


(18) 
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where  X^  =  2 nc/u>  is  the  wavelength  of  the  carrier.  The  expression  in  (18) 
is  based  on  the  assumption  that,  in  the  trapped  particle  regime,  electrons 
undergo  synchrotron  oscillations  in  the  ponderomotive  potential  wells. 

From  Fig.  7(a),  the  wavelength  of  the  sideband  radiation  is  within  10%  of 
that  given  by  Eq.  (18).  Coupling  of  a  parasitic  wavelength  A  ^  to  the 
carrier  at  Xg  via  the  synchrotron  oscillations  leads  to  a  maximum  growth  of 
the  instability  when  Eq.  (18)  is  satisfied.  In  the  presence  of  the 
sidebands,  the  optical  field  will  be  modulated  at  a  wavelength  XmQ^  given 

by 


1/2 


mod 


1  +  a‘/2' 

V 


2awlaolfB  ; 


(19) 


-4 

Taking  |aQ  |  «  0.8x10  ,  the  modulation  period  predicted  by  (19)  has  the 

approximate  value  0.12  mm,  which  is  in  close  agreement  with  the  period  of 
the  prominent  oscillations  observed  in  Figs.  4(a)  and  5(a).  The  close 
agreement  is,  however,  fortuitous  due  to  several  reasons.  For  one  thing, 
Eq.  (19)  is  derived  in  a  one-dimensional  calculation  for  the  case  where 
|aQ  |  is  constant  in  z,  whereas  there  is  significant  modulation  (~  100%)  of 
the  amplitude  in  Fig.  4(a).  Additionally,  for  the  finite  pulse  case 
examined  here  there  is  a  spread  in  the  wavelength  of  the  main  signal  and 
this  leads  to  the  rather  complicated  spectral  distribution  observed  in  Fig. 


7(a). 

The  asymmetric  profile  of  the  optical  pulse  in  Figs.  4(b)  and  5(b)  may 
be  understood  in  the  following  way.  At  the  wiggler  entrance,  the  electron 
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pulse  and  the  optical  pulse  are  symmetric  in  z  and  superimposed  on  top  of 
each  other.  Along  the  viggler,  the  electron  pulse  slips  behind  relative  to 
the  optical  pulse,  thus  tending  to  amplify  the  trailing  side  of  the  latter 
more  than  the  leading  edge.  Still  farther  behind  along  the  trailing  side 
of  the  optical  pulse,  the  field  amplitude  is  so  small  that  a  substantial 
fraction  of  electrons  are  not  trapped  in  the  ponderomotive  bucket.  This, 
in  conjunction  with  the  fact  that  diffraction  further  reduces  the  field 
amplitude  for  0  <  z  <  1  mm,  accounts  for  the  very  small  amplitude  in  this 
region. 

In  Fig.  8  (a)  and  (b)  the  extraction  efficiency  is  shown  for  the  two 
tapering  rates,  Cases  (a)  and  (b).  As  usual,  the  efficiency,  for  the 
slower  tapering  rate,  is  modulated  by  the  synchrotron  oscillations  of  the 
electrons.  For  the  more  rapid  rate  of  tapering  the  synchrotron 
oscillations  are  barely  noticeable.  Note  also  that  there  is  a  ten-fold 
improvement  in  the  efficiency  compared  to  the  slower  tapering  rate. 

-  Tapering  reduces  sideband  modulation  by  decreasing  the  trapping 

4  12 

fraction  and  by  distorting  the  synchrotron  motion.  ’  The  trapping 
fraction  drops  from  -  40X  in  Case  (a)  to  -  35Z  in  Case  (b).  A  measure  of 

4 

the  distortion  of  electron  orbits  is  given  by 


c(drr/dz) 


<2  (y-y  )> 
synv  ' r ' 


where 


(20) 


syn 


ck 

w 


(21) 


is  the  synchrotron  frequency.  Q  is  the  ratio  of  the  change  in  energy 
cdYr/dz  due  to  tapering  and  the  change  in  energy  <2(y  -  Yr)>  due  to 
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synchrotron  motion.  For  Case  (a),  Q  =  1%,  indicating  a  slight  distortion, 
whereas  for  Case  (b),  Q  =  25%,  indicating  significant  modification  of  the 
synchrotron  motion  and  thus,  reduced  sideband  modulation  as  is  indeed 
observed  in  Fig.  4(b)  and  5(b). 

(b)  Optical  Guiding 

Optical  guiding  in  the  context  of  FELs  has  been  the  subject  of 

numerous  discussions.*  Briefly,  the  resonant  interaction  between  the 

electrons  and  the  optical  field  tends  to  guide  the  radiation  along  the 

direction  of  the  electron  beam. 

3  10  19-25 

Two  causes  of  guiding,  ’  *  gain  focusing  and  refractive  guiding, 

have  been  distinguished  based  on  the  notion  of  a  complex  refractive 
24 

index.  In  general  these  two  participate  simultaneously  and  their 
combined  effect  on  the  spot  size  can  only  be  ascertained  via  the  envelope 

3 

equation  for  the  radiation  beam  in  an  FEL.  Refractive  guiding,  which  is 

described  by  the  reactive  (real)  part  of  the  refractive  index,  is  due  to 

10  20  21 

the  phase  shift  of  light,  as  implied  by  the  fact  that  vphase  <  c.  ’  ’ 

A  distinguishing  characteristic  of  this  type  of  guiding  is  that  the  optical 
wavefronts  are  plane  under  the  conditions  of  perfect  guiding.  The  other 
type  of  guiding,  gain  focusing,  is  described  by  the  resistive  (imaginary) 
part  of  the  refractive  index.  In  this  case,  under  the  conditions  of 
perfect  guiding, the  wavefronts  are  convex,  corresponding  to  the  fact  that 
there  is  a  net  power  flow,  due  to  diffraction,  away  from  the  electron  beam. 

Figures  9  (a)  and  (b)  show  the  spot  size  rg(z)  (dashed  line),  and 
wave-front  curvature  R~*(z)  (solid  line)  for  Cases  (a)  and  (b).  In  Fig.  9 
(a),  which  correspond  to  the  slow  taper,  one  discerns  the  modulations 
characteristic  of  sidebands,  as  in  Figs.  4  (a)  and  5  (a).  Besides  the 
modulations,  however,  Fig.  9  shows  a  dramatic  example  of  optical  guiding. 
Referring  to  Figs.  4  (a)  and  (b),  it  is  apparent  that  in  the  region  where 
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•■he  optical  field  amplitude  is  significant,  the  spot  size  of  the  radiation 
is  considerably  smaller  than  in  neighboring  regions.  In  fact,  in  Fig.  9  in 
the  regions  z  <  1.0  mm  or  z  >  3.0  mm  the  spot  size  rg  and  curvature  R_* 
take  on  values  corresponding  to  vacuum  diffraction  from  the  initial  values 
of  rg  and  at  the  viggler  entrance.  This  is  obviously  as  it  should  be, 
since  there  is  hardly  any  radiation  in  these  two  regions  (z  <  1.0  mm  and  z 
>  3.0  mm)  during  the  simulation. 

Another  significant  feature  of  Fig.  9  is  revealed  by  noting  that  the 
curvature  of  the  wavefronts  R-*  is  less  for  the  less  rapidly  tapered  case. 
In  particular,  in  the  vicinity  of  z  =■  1.2  mm  in  Fig.  9(a),  the  curvature  is 
negative,  indicating  that  the  wavefronts  in  this  region  are  in  fact 
concave.  On  the  other  hand,  in  the  case  of  rapid  tapering,  Fig.  9(b),  the 
wavefronts  are  convex  all  through  the  pulse,  indicating  a  flow  of  power 
away  from  electron  beam  throughout  the  pulse. 

To  shed  more  light  on  this  behavior,  it  is  useful  to  compare  the 
refractive  index  of  the  FEL  for  these  two  cases.  Using  a  refractive  index 


y  (r,z,t)  to  represent  the  entire  FEL  interaction,  the  vector  potential  may 


be  written  as 


3,10,20,21,27 


a(r,z, t) 


(22a) 


where 


(22b) 


p  *  1  +  (to^/w)  (ay/2|a|)  <exp(-i £.)/r> i 


(22c) 


is  the  electron  beam  plasma  frequency  and  <•••>  denotes  an  average  over 
the  local  particle  distribution.  In  terms  of  Fn  defined  by  Eq.  (2d), 
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.  t,  .  1  c  1,0  FnLn  ^2/rs2) 

"  ’Z  2  “  f  n  a  L  (2  ’  ' 

n=0  n  nl, 


|2r2/r_2] 


Figures  10  (a)  and  (b)  sb:>w  the  real  part  (Real(p)  -  1)  and  Figs.  11  (a) 
and  (b)  show  the  imaginary  part  (Imag.(p))  of  the  refractive  index  on  axis, 
r  =  0,  at  the  end  of  the  viggler  for  Cases  (a)  and  (b).  Comparing  Figs.  10 
(a)  and  (b)  it  is  apparent  that  there  is  more  refractive  guiding  in  the 
former  case,  which  is  the  less  rapidly  tapered  example.  This  is  consistent 
with  the  fact  that  for  slow  tapering,  £r  <<  1, 

cos  £ 

Rew  _  l  -  (23a) 

dominates  over  the  gain  term 
sin  l 

Imp  -  -  <  r>.  (23b) 


Equation  (23)  is  obtained  by  making  use  of  Eqs.  (2d),  (15)  and  (16).  On 
the  other  hand,  the  more  rapidly  tapered  case  of  Fig.  11  (b)  is  seen  to 
have  a  larger  net  gain  in  the  region  where  the  optical  field  is 
significant.  [Note  from  Eqs.  (22)  that  Imp  <  0  corresponds  to  gain.]  In 
the  same  region,  Fig.  11  (a)  indicates  that  the  net  gain  is  approximately 
zero  after  averaging  over  the  synchrotron  modulations.  Again,  this  feature 
of  Figs.  11  (a)  and  (b)  is  consistent  with  the  fact  that  in  the  latter  the 
tapering  rate  is  faster  and  therefore  the  gain  should  be  larger.  Moreover, 
it  is  seen  in  Figs.  9  (a)  and  (b)  that  the  larger  gain  corresponds  to  a 
more  convex  optical  wavefront. 

The  implication  of  these  results  with  regard  to  the  spot  size  may  be 

2  .3 

ascertained  by  a  consideration  of  the  terms  in  K  in  the  envelope  equation 
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rs'  +  K2(*,t.rs»laol)rs  =  °’ 

s  (2c/«)^(-l  +  2Ccos^r  +  C^sin^i^  +  (co/2c)r^C'sin£r)r~^, 

where  is  the  amplitude  of  the  fundamental  Gaussian  mode, 

C=(2I^/17rr)Hav/ |aQ  | ,  '  e  3/32  t-  c_*3/3t,  1^  is  the  beam  current  in  kilo- 

* 

Amperes,  H,  which  is  a  form-factor  related  to  the  transverse  profile  of  the 

electron  beam,  is  roughly  a  constant  and  close  to  unity  herein,  rr  is  the 

relativistic  factor  for  a  resonant  electron,  and  is  the  resonance  phase 

2 

approximation  for  <£> .  The  -1  in  the  expression  for  K  is  due  to  vacuum 
diffraction,  2Ccos£r  contributes  to  refractive  guiding  arising  from  the 
real  part  of  u,  the  third  and  the  fourth  terms,  due  to  the  imaginary  part 
of  y,  contribute  to  gain  focusing.  Taking  account  of  the  fraction  of 

trapped  electrons,  one  finds  that  in  going  from  Case  (a)  to  Case  (b)  the 

22  -3-2 

gain-focusing  tetir  C  sin  increases  from  2.4x10  to  7x10  .  The  other 

2  -3  -2 

term,  («/2c)rsC'sin£r ,  changes  from  -4.4x10  to  -2.4x10  ,  the  negative 

sign  indicating  a  defocusing  contribution.  On  the  other  hand,  the 

refractive  guiding  term  2Ccos£^_  decreases  from  3.1  to  1.6.  The  increase  in 

the  magnitude  of  the  gain  focusing  terms  in  going  from  Case  (a)  to  Case  (b) 

is  principally  due  to  the  increase  in  Concurrently,  the  50X  reduction 

in  the  refractive  guiding  term  is  due  to  the  increase  in  |ag|  and  the 

decrease  in  aw-  Since  the  refractive  guiding  term  is  the  dominant  term, 

2 

the  reduction  in  its  value  leads  to  a  decrease  in  K  ,  and  hence  to  reduced 
optical  guiding.  The  net  effect  is  the  increase  in  the  spot  size  and  the 
curvature  observed  in  Pig.  9(b)  as  compared  to  Fig.  9(a).  In  other  words, 
the  wavefronts  become  increasingly  convex  with  faster  tapering  rate. 

It  should  be  remarked  that  since  the  peak  output  field  amplitudes  for 
the  two  different  tapering  rates  are  about  equal  [cf.  Figs.  4  (a)  and  ( b ) ] , 
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the  increased  peak  output  power  of  the  more  rapidly  tapered  case  is 
principally  due  to  the  larger  spot  size.  Thus,  by  increasing  the  tapering 
rate  there  results  a  reduction  in  the  refractive  guiding;  however,  the 
increased  gain  then  causes  the  wavefronts  to  become  more  convex, 
diffracting  the  optical  field  into  a  larger  cross-section.  As  a  result  the 
output  power,  which  is  proportional  to  (|a|r  )  ,  increases  as  shown  in  Fig. 
3(b).  Detailed  examination  of  the  results  for  all  the  tapering  rates 
indicates  that  the  optical  field  amplitude  is  very  little  affected  by 
tapering  -  aside  from  reducing  synchrotron  modulations.  As  a  result, 
increased  output  power  and  extraction  efficiency  is  principally  due  to  the 
larger  lateral  extent  of  the  optical  field.  This  is  in  sharp  contrast  to 
the  situation  in  one  dimension  wherein  rg  is  necessarily  constant  and 
increased  efficiency  can  only  come  about  as  a  result  of  an  increase  in  the 
field  amplitude. 

Figure  12  summarizes  the  results  for  the  nine  tapering  rates  -dy^dz  = 
0.1^  0.3,  ...»  1.7  m-*,  corresponding  to  final  resonant  angles  =  1.8°, 
3.5°,  ...,  35°.  Beyond  -dy^dz  =  0.3  m""\  the  amplitude  |a|  is  fairly 
constant  up  to  -dy^/dz  «=  1.3  m”\  after  vhich  it  decreases.  However,  there 
is  a  near-monotonic  increase  in  the  spot  size.  Therefore,  it  is  the 
increased  transverse  extent  of  the  optical  field  —  and  not  an  increase  in 

intensity  —  that  is  responsible  for  the  enhancement  in  the  power  (« 

2 

|rsa|  )  observed  in  Figs.  3  (b)  and  5  (b).  Based  on  the  desired  output 
power  and  the  constraint  on  the  maximum  spot  size  one  can  determine  the 
optimal  tapering  from  Fig.  12. 
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V. 


Conclusions 


This  paper  presents  the  results  from  the  simulation  of  a  high-power, 
high-current,  short  wavelength  FEL.  The  simulations  include  many  of  the 
important  aspects  of  a  realistic  experiment,  including  finite-pulse 
effects,  diffraction  and  emittance.  A  wide  range  of  wiggler-f ield  tapering 
rates  have  been  examined.  Summarizing  the  results,  it  is  found  that  in 
addition  to  enhancing  the  efficiency,  tapering  improves  the  quality  of  the 
output  by  suppressing  sideband  growth.  Further,  it  is  found  that  as  the 
tapering  rate  is  increased  there  is  a  gradual  transition  from  a  refractive- 
guiding  regime  to  one  where  gain  focussing  dominates,  with  optical  power 
diffracting  laterally  along  the  convex  wavefronts.  The  increased  transverse 
extent  of  the  optical  field,  rather  than  an  increase  in  the  field  amplitude 
is  the  major  reason  for  efficiency  enhancement  as  the  tapering  rate  is 
increased. 
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Table  I.  Parameters  for  a  high-power,  rf-linac  FEL 


Electron  Beam 

2 

Energy,  ymc 
Current,  1^ 

Normalized  edge  emittance,  ex  =  ey 
Radius,  rbo 

Betatron  period,  2Ji/k^ 

Pulse  length 
Viggler 

Induction,  Bv 
Period,  2rt/kw 
Length 

Input  Radiation 

Wavelength,  2nc/co 
Spot  size,  rg(0 
Rayleigh  range,  7^ 

Radius  of  curvature  of  optical  wavefronts,  Zg/a(0) 
Pulse  length  (FWHM) 

Peak  input  power 


175  MeV 
50  A 

153  mm-mrad 
1  mm 
11.5  m 

6.7  ps 

6.4  kG 

4.7  cm 
42  m 

1  ym 
1.25  mm 
4  m 
10  m 

31.4  ps 
450  MV 
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Table  II.  Resonance  phases  and  tapering  rates  (The  energy  spread  due  to 
emittance  is  approximately  the  same  for  all  tapering  rates) 


dy  /dz 
r 

(m-1) 

Initial 

Resonance  Phase 

lr  (°) 

Full  bucket 

height 

[  ^rr'/rrl  bucket 
(Z) 

-0.1 

2.6 

0.96 

-0.3 

7.9 

0.89 

-0.5 

13.2 

0.82 

-0.7 

18.7 

0.75 

-0.9 

24.3 

0.68 

-1.1 

S  1.2 

0.6 

-1.3 

36.5 

0.5 

-1.5 

4!. 2 

0.44 

-1.7 

48.0 

0.35 

Energy  spread 
due  to 
emittance 

1  Wendt 
(Z) _ 


0.12 
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DISTANCE  ALONG  WIGGLER  (m) 

Tapering  of  normalized  viggler  amplitude  along  the  wiggler. 
Radiation  Power  in  GW  at  viggler  exit,  (a)  slow  tapering  rate 
dyr/dz  =  -0.1  Case  (a);  (b)  rapid  tapering  rate,  drr/dz  = 

-1.3  oT*,  Case  (b). 
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Fig.  4  Vector  potential  |a(z,r«0)|  and  linear  electron  density  N(z)  at 
viggler  exit,  (a)  slov  taper,  Case  (a);  (b)  rapid  taper,  Case 

n 

(b).  Note  that  intensity  «  |a|  . 
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55.0 


Fig.  5 


Radiation  Power  in  GW  at  wiggler  exit,  (a)  slow  tap 
(b)  rapid  taper,  Case  (b). 
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.  8 


DISTANCE  AL 

Extraction  efficiency  alon 
(a);  (b)  rapid  taper,  Case 


NORMALIZED 

CURVATURE  c,  SPOT  SIZE  r  (cm)  AMPLITUDE  |a|  PEAK  POWER  (GW)  EFFICIENCY  (%) 


-dyrldz  (m_1) 
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